Effect of Spinless Impurities on Reduction of T c in High T c Superconductors 
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The notion of a finite pairing interaction energy range suggested by Nam, results in some states at 
the Fermi level not participating in pairings when there are scattering centers such as impurities. 
The fact that not all states at the Fermi level participate in pairing is shown to suppress T c in an 
isotropic superconductor and destroy superconductivity. We have presented quantitative calcula- 
tions of T c reduced via spinless impurities, in good agreements with data of Zn-doped YBCO and 
LSCO, respectively. It is not necessary to have the anisotropic order parameter, to account for the 
destruction of superconductivity via non-magnetic impurities. 
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One of intriguing experiments in high T c superconduc- 
tors (HTS) is that the non-magnetic impurity Zn reduces 
T c and destroy superconductivity in Zn-doped YBCO 
H I and Zn-doped LSCO @, respectively. 

The spinless impurity scattering, according to Ander- 
son Q , would not change thermodynamical properties of 
superconductor, such as T c , providing the order param- 
eter being a uniform function, that is, a constant. He 
argued that the pairing of the state and its time reversal 
state in the scattering quantum space, would yield the 
same order parameter as that in the unscattered case, 
since the states in the former space can be mapped via 
the unitary transformation from those in the latter. In 
practice, one of crucial parameters to determine T c is the 
BCS coupling parameter 



g = (N(0)Vi 



BCS/ 



(1) 



with the density of states at the Fermi level N(0) and 
the BCS pairing interaction Vbcs averaged. We expect, 
thus, that the anisotropic pairing interaction would af- 
fect T c . In fact, Markowitz and Kadanoff || have ex- 
amined the anisotropic effect on T c and obtained good 
agreements between calculations and data for low con- 
centrations of non-magnetic impurities in low T c super- 
conductors (LTS). 

On the other hand, Abrikosov and Gor'kov (AG) 0] 
have shown that the spin-flip scatterings suppress T c and 
destroy superconductivity. Also, when the order parame- 
ter has nodes, say p- or d-wave pairing states, Hirschfeld, 
Wolfle, and Einzel M argued that the resonance impurity 
scattering causes T c be reduced in the way of AG formula. 
In essence, the anisotropic effects and pair breaking scat- 
terings, such as, spin-flip scattering, are known to be, up 
to now, the causes for destruction of superconductivity. 
We will show shortly a new chapter of destruction of su- 
perconductivity via non-magnetic impurity scatterings in 
an isotropic superconductor. 

The notion of a finite pairing interaction energy range 
Td H results in the incomplete condensation in which not 
all states participate in pairings. The states not partici- 
pating in pairings yield low energy states responsible for 



the linear T dependence of superelectron density at low 
T in an isotropic superconductor ||l0| . The quantitative 
calculations of magnetic penetration depth length in all 
T ranges based on a finite Td are in good agreements 
with data of YBCO jjjj , BSCCO ||, HBCCO ||, 
LSCO @, and Sr214 respectively. Moreover, the 
incomplete condensation results in the picture of mul- 
ticonnected superconductors |l8| ] which can account for 
the 7r-phase shift in the Pb-YBCO SQUID || and \ 
magnetic fluxoid quantum observed in the YBCO ring 
with three grain boundary junctions pQ] . 

First we recapitulate the pertinent results of a finite Td 
To see the phase transition, T c should have a finite 
value, that is, neither zero nor infinite. To have a finite 
value T c , Td should be finite, since T c is scaled with Td 
within the pairing theory. In other words, we may write 
the order parameter A(k,tv) as M 



A 





for 
for 



<T d 
>T d 



(2) 



for all dynamical energies uj. Here tk is the usual normal 
state excitation energy of wave number fc, measured with 
respect to the Fermi energy. Hereafter units of % = c = 
ks = 1 are used. And later A(T) for A may be used as 
well. The A is a solution of BCS like equation 



de ( E 
— tanh — — 
E \2T 



with E = (e 2 + A 2 ) 1 / 2 . For A(T c0 ) = 0, we get 
l/. 9 =(2/ 7 r)^(2/j)tan- 1 ( 2/ /j), 



(3) 



(4) 



where y = Td/{^T c o) with the transition temperature 
T c q for a pure system, and sum is over the positive odd 
integers j. The factor of arctangent function makes the 
sum converge. For large y, Eq. (^) yields the BCS result 
Tco(BCS). The quantitative calculations of T c q and A 
vs g and T are given in Ref . 1 1 1 . The T c q from Eq. 
(0) is always greater than T c0 (BCS) and does not have 
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any upper limit. The fact is that for large g > 2.32, 
T c o increases with increasing g as T c q ks gTj/2. One 
interesting value of g = 0.657 yields T c q ~ 100 K with 
Td = 400 K which is of the order of Debye temperature 
in HTS. This value of g may be realized in YBCO by 
considering the electron-phonon interaction of the order 
of X p = 1.3 - 2.3 jn|. 

The key idea of this letter is a novel one. When 
there are scattering centers such as impurities in a Fermi 
system, the spectral weight distribution function at the 
Fermi level, that is, the imaginary part of the Green's 
function at u> = 0, may be considered as a Lorentzian 
form. The sum of spectral weights outside Td < \ek\, 
that is, the integral of the imaginary part of the Green's 
function with t\. in the ranges of |ejt| > Td, is given as 



N{Q)R = iV(0)(2/ 7 r)tan- 1 (r/T ti ), 



(5) 



where T is the imaginary part of self-energy due to the 
impurity scatterings, or the half impurity scattering rate 
l/(2r). The states of N(0)R do not participate in pair- 
ings, and are predicted [|22| to be responsible for the lin- 
ear T dependence of the specific heat at low T in Zn- 
doped YBCO and also for the T 2 term in the 

magnetic penetration depth length at low T in the im- 
perfect HTS such as Sr214 @- The effective density 
of states at the Fermi level participating in pairings, be- 
comes AT(0)(1 — R) less than the usual one N(0), and re- 
sults in the reduced T c . For simplicity, we assumed here 
the effective BCS coupling parameter g be not changed 
via impurity scatterings. Also the Born approximation 
is used for calculations of impurity scatterings. 

To incorporate the above key idea into the calculation 
of T c , we use the Green's function formalism p5j . In 
the thermal Green's function scheme, i.e., replacing uj 
by iu> n = i(2n + 1)ttT with integer n, we may, taking 
into account the order parameter of Eq. (|^) consistently, 
rewrite Eqs. (3.5a, b, c) of Ref. |^5| as 



Z n u> n = uj n + (r + r s )[$i(w n ) + $ 2 (w„)], 

z n A n = A + (r-r s )$ 3 K), 

71 

*i(w„) = (u Jn /E n )(2/ir)ta,n- 1 [T d /{Z n E n )}, 
$ 2 K) = (2/n)ta3ar 1 {Z n u n /T d ), 
$3(w„) = (A n /u) n )$i(u n ), 



(6) 
(7) 
(8) 



Here T and T s are the imaginary parts of self-energy via 
non-magnetic and magnetic impurity scatterings, respec- 
tively. When r = = T a , Z n = 1, and A„ = A, Eq. (§) 
becomes Eq. (||) of the BCS like equation, as it should. 

For T -> T c , A -> 0, from Eqs. (§), (0), and (|h, we, 
after simple algebras, get the equation for T c as p2fl 



1 



E2irT c r n 



(9) 



D„ 



r n = (2/tt) tan-^lVK + T + T s )], 



u> n = (2n + 1)ttT c . 

The sum is over integers n > 0. For T = = T s , Eq. 
(0) becomes Eq. ([|) for a pure superconductor. In the 
infinite Td limit, equivalently small g limit, that is, LTS 
case, combining Eq. (^) and Eq. (|h yields the AG result, 
with r n = 1 and D n — 2T S , as |Q,g5| 



ln(T c /T c0 ) = 



2n+ 1 



2n+ 1 



(10) 



where the depairing parameter d = 2T s /(ttT c ). Further- 
more, T c becomes independent of T, and the Anderson 
theorem Q follows with T s = 0. In general cases, com- 
bining Eq. (|^) and Eq. (^|), we calculate easily T c /T c o 
and T c /Td vs T, T s , and g, respectively. The calculations 
for g = 0.3, 0.45, 0.657, 1, and 2, have been carried out, 
and some results are presented here. 

In Fig. 1 is shown T c /T d vs 2r/T d and 2T S /T d for 
g = 0.657 for example. The general shape of T c /Td 
appears to be similar in any direction of T and T s . 




2r/Td a6 



2rVTd 

FIG. 1. T c /T d vs 2F/T d and 2T S /T d for g = 0.657. 
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To see somewhat in detail, in Fig. 2 are shown T c /T d 
vs 2Y/T d for T s = and T c /T d vs 2T S /T d for T = 
0, for g = 0.3, 0.45, 0.657, 1, and 2, respectively. 
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FIG. 2. The left panel is T c /T d vs 2T /T d for T s = 0, while 
the right panel is T c /T d vs 2r s /T d for r = 0. 

The reductions of T c via Y are slow for small g, and the 
steepness increases with increasing g. On the other hand, 
the reductions of T c via T s are almost the same shape for 
all g values. 

To see some more differences between non-magnetic 
and magnetic impurity scatterings, we have changed the 
scales of 2r/T d and 2T S /T d by 2r/A and 2r s /A , 
with Ao = A(0) the order parameter in a pure super- 
conductor, at T = 0, and a function of g. In Fig. 
3 are shown T c /T d vs 2r/A and T c /T d vs 2r s /A . 




FIG. 3. The left panel and the right panel are T c /T d 
vs 2r/A with T s = and T c /T d vs 2r s /A with T = 0, 
respectively. The normalized T c /T c o are given in the insets. 



The insets in Fig. 3 contain the normalized T c /T c o vs 
2r/A and T c /T c0 vs 2r s /A , respectively. The normal- 
ized T c /T c0 for all cases have almost the same shape. 
In the T c = limit, we may rewrite Eq. (||) as ]22| 

1 f°° , rJx) , , 

g J x + D c {xY 

r c (x) = (2/tt) tan- 1 [{T d /A )/(x + x 1 +x 2 )], 

D c (x) = [1 - r c (x)]xi + [1 + r c (x)]x2, 

an = r c /A , 

x 2 = r sc /A . 

The phase diagram for T c = is shown in Fig. 4. 
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FIG. 4. The phase diagram on V — F s plane. The critical 
values of F c and F sc are given in the inset. 

In the inset in Fig. 4, the critical vaules of 2r c /Ao de- 
crease with increasing g, contrast to those of 2T C /T d in 
Fig. 2 increase with increasing g 7 since T c increases with 
increasing g. However, both 2r sc /Ao and 2T SC /T d in- 
crease with increasing g. For small g limit, LTS case, 
the AG result 2r sc /Ao = 0.5 is obtained. In the large g 
limit, HTS case, we get from Eq. ([□]) as f| 

2r c /A + 2r sc /A - 4/tt. (12) 

This can be realized from Eq. (|^) as well by observ- 
ing r n — and D n — T + T s for large g or small T d 
limit. This implies the effects of non-magnetic and mag- 
netic scatterings on T c are the same in HTS. Data in 
HTS indeed indicated such as shown in Fig. 5. To com- 
pare calculations with data, we used T c /T c q vs T/T c in 
Fig. 5. In addition to data of Zn-doped HTS, those of 



3 



Ni-doped YBCO and LSCO also are included to demon- 
strate both non-magnetic and magnetic scatterings af- 
fect T c in the same way for large g, that is, HTS case. 
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YBCO(Zn) Ref [1] 
YBCO(Zn) Ref. [2] 
YBCO(Zn) Ref. [3] 
YBCO(Ni) Ref. [3] 
LSCO(Ni) Ref. [3] 
LSCO(Zn) Ref. [3] 
LSCO(Zn) Ref. [4] 
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g=0.45 
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FIG. 5. Comparison between calculations from Eq. (£j) and 
data of Zn-doped and Ni-doped YBCO and LSCO. Only two 
theoretical curves (g = 2.0 and 0.45) are shown. Others are 
bounded by these two curves. 



The data of YBCO || are poor but that of YBCO 
are in good agreements with our calculations. The data 
of LSCO § are satisfactory. In all, data are understood 
in the framework of a finite T<j. 

In summary, even though the model of the order pa- 
rameter of Eq. (^|) is ideal, and the Born approximation 
of impurity scatterings in an isotropic superconductor is 
used, the calculations of T c vs non-magnetic and mag- 
netic impurity scatterings, are in good agreements with 
data of YBCO and LSCO, respectively. The notable re- 
sults are the normalized T c /T c q vs F/r c are almost the 
same each other for all g values, and T c /T c q vs r s /r sc 
yields the AG result. For improvement, one may use 
the T matrix formulation for impurity scatterings, and 
take into account the retardation effect on the order pa- 
rameter. However, we expect the T c /T c q would not be 
changed much from that discussed here. The spinlcss 
impurity scatterings can destroy superconductivity in an 
isotropic superconductor as well. The notion of a finite 
Td is sound theoretically. We hope it would account for 
other unsolved problems in HTS. 
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